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Abstract 

We present an analysis of the decay K~^7r~ based on FOCUS data. From a sam- 
ple of 234 signal events, we find a branching ratio of ppoI|^^-^+| = (0-429 Iq qq^ =b 
0.027)% under the assumptions of no mixing and no CP violation. Allowing for CP 
violation, we find a branching ratio of (0.429 it 0.063 it 0.028)% and a CP asym- 
metry of 0.18 lb 0.14 lb 0.04. The branching ratio for the case of mixing with no CP 
violation is (0.381 lQ^g3 ib 0.092)%. We also present limits on charm mixing. 

Key words: 

PACS: 13.25.Ft 14.40.Lb 12.15.Ff 



1 Introduction 

While mixing in the strange and beauty quark sectors has been observed 
for many years, charm mixing remains elusive. In the Standard Model the 
charm mixing rate is greatly suppressed by small CKM matrix elements and 
strong GIM suppression. It is this very fact, however, that provides a unique 
opportunity to search for new physics. 
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This paper represents a continuation of mixing studies by FOCUS begun 
with a measurement of ycp looking for a hfetime difference between CP 
even and mixed CP states [1] and continued with an analysis of the decay 
— > K'^7r~ [2]. A summary of charm mixing can be found in the recent 
Particle Data Group review [3] . 

The search for charm mixing requires tagging the flavor of a neutral charm 
meson at production and tagging the flavor again at decay The production 
flavor is determined using D*'^ — > D^nf decays. The charge of the slow pion 
(tTs) determines the flavor of the produced neutral D meson. The flavor of the 
decaying is determined by reconstructing a Cabibbo favored (CF) decay 
with a charged kaon such as K~ti^ . A right-sign (RS) decay is defined 

as one in which the kaon and soft pion charges are opposite while same kaon 
and TTs charges are wrong-sign (WS) decays. The can also decay directly to 
K'^T^~ via a doubly Cabibbo suppressed decay (DCSD). In the small mixing 
limit (known to be true experimentally), the time dependent wrong-sign decay 
rate can be written as: 

i?ws(t) = e-^* [Rd + ^DV'Tt + ^ (x" + y'^) VH'^^ (1) 

where the three terms correspond to DCSD, interference between mixing and 
DCSD, and mixing. Rd is the DCS branching ratio relative to the Cabibbo 
favored mode. The mixing parameters x' = xcos6ktt + ysindxTT and y' = 
y cos Sktv — X sin 6k-k are rotated versions of the mixing parameters x = 
and y = ^ which are the mass and lifetime splitting terms, respectively. 
The angle Sktt is the strong phase between the CF and DCS decay, generally 
expected to be small. Analogously to the DCSD rate, Rm = (x''^ + y''^^ /2 
is the time independent mixing rate. Discovery of hadronic charm mixing 
requires separating the three components of Eq. 1 using the different lifetime 
distributions. Note also that Eq. 1 only has a x'^ and therefore the sign of x' 
cannot be determined from this analysis and the relevant fit variable is x'"^. 

CP violation can in principle occur in any of the three wrong-sign components 
resulting in three additional terms, Ajj, Am, and (j). A^ and Am are CP 
asymmetries associated with the DCSD and mixing terms, respectively. The 
CP violation term associated with the interference simply rotates Skw by ±0. 
One choice for the CP violating time dependent wrong-sign rate is: 





[y cos ^ X sin 0) Ft (2 j 
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where in all cases the upper (lower) sign refers to an initial D° ■ We can 

use Eq. 2 on the full D° + data set to directly extract the six parameters 
of interest, Ro, x''^, y', Ajj, Am, and 0. There is a four-fold degeneracy in 
Eq. 2. A 2-fold degeneracy appears when x\ y', and (f) are replaced with —x', 
—y', and (p + n which does not change Eq. 2. Requiring |0| < 7i/2 removes 
this degeneracy and follows the convention of BABAR [4] . Operationally, this 

is done by replacing cos0 with \/l — sin^ cf) and fitting for sin0 rather than 
(f). The other 2-fold degeneracy occurs when replacing x' and (f) with —x' and 
— 0, again leaving Eq. 2 unchanged. We remove this degeneracy by replacing 
x' with \x'\. This convention avoids a specious dependence on the sign of x'. 

An alternative method for extracting the six parameters of interest is to use 

Eq. 1 on the and D° samples separately and obtain values of Rd, x''^, 
and y' for each. We obtain |i?J,a;'^ '^'^} ^^'^ {-^D)^'^ } ^^'^^ 
and samples, respectively. With = ^x'^ -|- y"^ ') /2, the CP violation 
asymmetries can be found as: 

-R^ - Rp , . _ Rm - Rm /on 

while a;'^ and y"^ are related by: 
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(5) 



Choosing either the -|- or — equations from Eqs. 4,5, we can square both sides 
and add the equations together. The result, along with rewriting Eq. 5, gives 
three relations: 
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Again replacing cos with y 1 — sin^ 0, one can solve these equations for x''^, y', 
and sin 0. There are two solutions corresponding to the relative sign between 



x''^ and x' . In this framework it is natural to define Rd = y RdRd- 
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2 Event reconstruction and selection 



The FOCUS experiment recorded data during the 1996-7 fixed-target run 
at Fermilab. A photon beam obtained from bremsstrahlung of 300 GeV elec- 
trons and positrons impinged on a set of BeO targets. Four sets of sihcon 
strip detectors, each with three views, were located downstream of the tar- 
gets for vertexing and track finding. For most of the run, two pairs of silicon 
strips were also interleaved with the target segments for more precise vertex- 
ing [6] . Charged particles were tracked and momentum analyzed as they passed 
through up to two dipole magnets and up to five sets of multiwire proportional 
chambers with four views each. Three multicell threshold Cerenkov counters, 
two electromagnetic calorimeters, and two muon detectors provided particle 
identification. A hadronic trigger passed 6 billion events for reconstruction. 
The average photon energy of reconstructed charm events is 175 GeV and the 
average momentum for this analysis is 75GeV/c. 

A candidate driven vertexing algorithm is used to reconstruct charm. In the 
case of K'n'^, two oppositely charged tracks are required to verticize 

with CL > 2%. The momentum and vertex location arc used as a "seed" 
track to find the production vertex which must have CL > 1%. The fiavor 
of the produced D meson is determined using the decay D*~^ D^^^t ■ The 
soft pion must be consistent with originating from the production vertex and 
the track is refit using the production vertex as an extra constraint. The en- 
ergy release {Q{D*) = M{D*) - M(D^) - m^) must be less than 55MeV/c^. 
Separating charm from hadronic background is primarily accomplished by re- 
quiring the decay vertex be distinct from the production vertex. A loose cut 
of £/(T£ > 2 is appUed where ^ is the distance between the two vertices and 
(7^ is the calculated uncertainty (<(j£> ~ 500 //m). Since hadronic reinter- 
actions can fake a decay, requiring the secondary vertex to be located out- 
side of target material reduces background. The out-of-material significance 
(Tout is positive (negative) for a vertex outside (inside) material. We require 
2£/cr£ + max (—2, (Tout) > 6. To ensure the decay tracks do not originate 
from the production vertex, a cut is made on the change in production vertex 
confidence level when either decay track is added to the vertex (ACLpri). 
We require ACLpn < 20%. Since the signal contains a charged kaon, infor- 
mation from the three Cerenkov counters effectively suppresses backgrounds. 
The Cerenkov identification algorithm [5] returns negative log-likelihood val- 
ues Wi(j) for track j and hypothesis i e {e,7r, In practice, differences 
in log- likelihoods between hypotheses are used such as AW^i^- = — W^-. 
We require A>V^A'(i^) > 0.5, ^WkM > -3, Wmin(7r) - W^(7r) > -5. and 
yVmhii'^s) - yVniT^s) > -5 where Wmin = min{Wi(z{e,7r,K,p})- We also require 
AW-^KiK) + AWkA^) > 3- To suppress doubly misidentified K-71+ 
decays (reconstructed as — > 7r~ir+), we require AWT^KiK) + A'Wktt{'^) > 
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8.5-0.5 [m^^(D0)1^ where M^^(D°) is the normahzed reflected mass obtained 
by switching the rest mass values for the tt and K, subtracting the nominal 

mass, and dividing by the calculated error. The remaining selection cri- 
teria are very efficient and mildly suppress some backgrounds. To suppress 
semileptonic decays we require the tt candidate not be identified as a muon 
by the muon detectors. Both pions must be identified as not being an electron 
by cither the Cerenkov system (A>Ve7r(7r) > 1) or a calorimeter. All tracks 
used in the analysis must be inconsistent with the beam direction to suppress 
tracks from 7^e"'"e^. Momentum requirements of p{K),p{7r) > 6GeV/c and 
p{ns) > 2GeV/c arc set. The calculated uncertainties on the lifetime and 
mass must be less than 150 fs and 25MeV/c^, respectively. The momentum 
asymmetry, ||f]5g|, must be less than ^^^jotceV/f^' • The summed p| of 

daughters with respect to the D° momentum vector must be greater than 
0.25 GeV^/c^. These last two cuts favor a decay of a heavy particle over combi- 
natoric background. If, after passing all these cuts, more than one D* candidate 
is found for one D^^ candidate, an additional cut may be made. If a right-sign 
(wrong-sign) candidate is found with 2.5MeV/c^ < Q{D*) < 9.5MeV/c^ 
then no wrong-sign (right-sign) candidates are allowed. 



3 Mixing fit 



A three dimensional binned maximum likelihood fit is used in this analysis. 
Two of the dimensions, M(D°) and Q{D*), are useful for separating signal 
from background while the third dimension, t{D^)^ is necessary to distinguish 
DCSD, mixing, and interference contributions. The mass is divided in 30 
equal bins from 1.75 to 2.05 GeV/c^. The energy release is divided into 110 
equal bins from to 55 MeV/c^. The proper hfetime is divided into 5 bins 
of 0.20 ps from to Ips, one bin from 1.0-1.25 ps, one bin from 1.25-1.75 ps, 
and one bin from 1.75-3 ps. These bins are much larger than our proper time 
resolution of ~35 fs. The right-sign and wrong-sign data are fit simultaneously. 
The negative log-likelihood £ = — log £ is simply the sum of the negative log- 
hkelihoods of each bin, (. = Yh""""" (-i- MiNUlT [7] is used to minimize £ and 
1-0" errors are defined as the point where I changes by 0.5 with respect to the 
minimum value using the MiNOS [7] algroithm. 

The fit has many components. These components generally have a 3-D shape 
determined by a Monte Carlo simulation and a yield which is free to fioat, some 
with weak constraints imposed. The components making up the fit which are 
modeled by the Monte Carlo are RS signal, WS signal (DCSD, interference, 
and mixing), real ^K~ti^ decay with a fake soft pion, reflections to both 
RS and WS {K-K+, 7r-7r+, 7r-7r+7r°, K^ti-ti+), RS background {K-(i+u and 
K~T^^T^^), and WS background (double misidentification of K~l'^v, X^tt+tt^, 
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K-Tr+). The 3-D Q{D*), & (r(L'°)) shapes are obtained from a 

Monte Carlo simulation with the PDG [3] value for the lifetime (appropri- 
ately modified for WS signal mixing and interference terms). The remaining 
background, labeled random, is some combination of non-charm and poorly 
reconstructed charm events. In the default fit, this contribution is modeled 
with functional forms in mass aexp(6m), energy release aq^^'^ + Pq^^^, and 
lifetime exp (— t/ri) +ricxp {—t/T2) where a, 6, a, /5, 77, ri, and T2 are free pa- 
rameters of the fit. In one fit variation, this contribution is modeled instead as 
a combination of minimum bias and generic charm events (without the decay 
modes accounted for above) from Monte Carlo. 

Penalty terms are added to the likelihood to ensure that many of the back- 
grounds described above are consistent with known branching ratios. The af- 
fected background components are the D° decays to K^K~, tt+tt", 7r+7r~7r°, 
ir°7r+7r~, K'T^^T^^, and K~l'^v. These six reflections appear in both right-sign 
and wrong-sign. Each of these 12 yields is a free parameter of the flt with a 
penalty term of 

[Ffit (L>°^ref ) - Fo (^°^ref )]' 

Penalty = 0.5 — 5 ^ 9 

[a(Fo(i^°^ref))f ^' 

added to the log-likelihood for each yield. Igt is the fitted yield, Yq is the 
expected yield based on the number of right-sign signal events, the relative 
efficiency between signal and the refiection, and the branching ratio between 
signal and the reflection. The uncertainty on the expected yield, cr(Fo), is 
based on the PDG [3] branching ratios inflated to account for effects such 
as Monte Carlo misidentification mismatch and efficiency variation due to 
mismodeled resonant substructure. The branching ratio and error used in the 
fit are shown in Table 1. The use of the penalty terms does not significantly 
affect the results as the signal shape is very different from these background 
shapes. The benefit is that correlations between backgrounds are significantly 
reduced resulting in a better behaving fit. The amount of background due 
to double misidentification of the right-sign signal mode is fixed based on 
the fitted right-sign yield and relative efficiency obtained from Monte Carlo. 
This background accounts for only 132 events in the entire data sample, 19 
in the signal region. The error on this number, determined from the Monte 
Carlo fidelity in estimating the K~K^ contribution, is less than 12%, 

resulting in a negligible contribution to the systematic error. 

During the course of this analysis it was observed that the Monte Carlo failed 
to reproduce a low mass tail from the signal >ir~7r+ decays. The program 
Photos [8] was used to model QED internal bremsstrahlung processes which 
occur during the charm decay. This addition produces a Monte Carlo signal 
shape which agrees with the data. Correlations in the Monte Carlo shapes 
between the three variables (M(Z)°), Q{D*), r(D°)) were examined. While 
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Table 1 

The PDG04 branching ratios and the values used in the fit to constrain the various 
background contributions. Errors are inflated to account for additional uncertainties 
such as particle identification and resonance substructure. For 7r~7r+7r° the errors 
are infiated due to a discounted measurement of (3.9 ih 1.1)%. 

Decay PDG BR (%) Used BR (%) 

r{K-K+)/r{K-7T+) 10.231^:27 10.26 ±0.32 

r(7r-7r+)/r(K-7r+) 3.62 ±0.10 3.58 ± 0.20 

r(7r-7r+7rO)/rtotai 1-1 ± 0.4 ± 0.2 1.1 ± 1.5 

r(K-7r+7r°) /T{K-7r+) 342 ± 22 342 ± 55 

r{K~e+iy) /Ftotai 6.76 ± 0.28 6.76 ± 0.90 

r(;^07r+7r-)/rtotai 5.97 ±0.35 5.97 ±0.875 

the lifetime was independent, the mass and energy release had significant 
correlations. Therefore, the Monte Carlo shapes are represented as a 2-D shape 
{M,Q) times a 1-D shape (r). The Monte Carlo shapes for individual modes 
are more than 50 times larger than the data sample. 



4 Systematic checks and uncertainties 



Mini Monte Carlo tests of the fit were used to verify the accuracy of the re- 
ported fit errors and to check for possible biases in the fit. 1000 independent 
data samples were generated from the best fit to our data. Each of these data 
samples were fit as real data. Distributions of fitted values indicate no signif- 
icant bias and accurate fit errors. Confidence levels based on the distribution 
of — logC and were found to be 63% and 27%, respectively. A direct mea- 
surement of the confidence level from the is impossible due to very few 
entries in most bins. 

Fit variants with different binning, different constraints, and different account- 
ing of the random background were tried. No significant differences were ob- 
served indicating a consistent fit with no unaccounted-for resolution effects. 
Variations of all the selection criteria were analyzed with no significant dif- 
ferences in the results. This is to be expected since the dependence on the 
Monte Carlo simulation is rather weak. For the branching ratios and asym- 
metries reported, systematic uncertainties due to fit variants and cut variants 
were obtained from the r.m.s of the variations and then added in quadrature. 
95% CL limits on the asymmetries were obtained by scaling the 1-a errors. 
When the asymmetry was measured directly in a fit, the scaling was deter- 
mined by the ratio of errors obtained from MiNOS for A log >C = 0.50 (l-cr) and 
AlogvC = 1.92 (95% CL). The scahng is calculated separately for positive and 
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negative errors. When the asymmetry was calculated from other parameters 
the naive scaling of 1.96 between l-cr and 95% CL uncertainties was used. 



5 Fit results 

We perform five types of fits to extract branching ratio, mixing, and CP vio- 
lation information. These fits are labeled A-E. Fits which constrain x'^ to lie 
in the physical region (> 0) are "primed." 

Fit A assumes no mixing and no CP violation (a:'^ = |/' = in Eq. 1). The 
wrong-sign branching ratio is measured to be -Rws = (0-429 ^q'qq'^ ± 0.027)%. 
Fit B assumes no mixing but allows for global CP violation (x'^ = y' = (J) = 
Am = in Eq. 2). The wrong-sign branching ratio and CP asymmetry arc 
found to be i?ws = (0.429 ± 0.063 ± 0.028)% and Ad = 0.18 ± 0.14 ± 0.041. 

Fits C and C allow for mixing without CP violation (Eq. 1). Representative 
plots of Fit C arc are shown in Figs. 1-3. Figure 1 show the right-sign projec- 
tions onto M(D°) and Q{D*), with and without a cut on Q{D*) and M{D°), 
respectively. Figure 2 shows the same projections for the wrong-sign events. 
Figure 3 shows the signal component of the wrong-sign lifetime distribution 
split into the three components, DCSD, mixing, and interference. 

The DCS branching ratio is found to be Rd = (0.381 ± 0.092)%. The 

negative log-IikcIihood is minimized at values of x'"^ = —0.059% and y' = 1.0% 
for Fit C. Due to significant correlations between x' and y' and the fact that 
< is unphysical, the interesting physics result is a contour in the x' , y' 
plane. The minimum negative log-likelihood with x''^ > occurs at x'^ = and 
y' — 0.5%. Although we fit for x'^, we choose to plot x'. The 95% CL contour 
is defined as the location in the x', y' plane where the change in log-likelihood 
reaches 3.00 relative to the minimum negative log-likelihood in the physical 
region of the x', y' plane. This method has been checked with a frequentist 
method using mini Monte Carlo with identical results. Systematic checks were 
performed with 120 fit and cut variants. The contour variation is consistent 
with differences in the returned value of x''^ and y' . To assess the systematic 
uncertainty, we first determine the change in negative log-likelihood between 
the global minimum and the x\ y' location for each of the 120 variants. We 
then find the value greater than 95% of these differences which is 0.48. We find 
the contour at which the change in negative log-likehhood is 3.00-1-0.48 = 3.48 
and call this the 95% CL including systematic uncertainties. The contours for 
Fit C are shown in Fig. 4. We define 95% CL limits on x' and y' based on 
the projection of the contour onto the respective axis. The 95% CL limit on 
Rm is defined as the maximum Rm value on the contour. The improvement in 
— \ogC compared to Fit A is 0.089 (0.082) for free (restricted) values of x''^ . 
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Fig. 1. Projection onto M(L>°) and Q{D*) of both the fit (Fit C) and data for 
right-sign events. The left plots are for all events while the events in the right plots 
have cuts in the non-plotted variable. Due to the very large signal-to-background, 
semi-log plots are shown. 

The last two fits allow for both mixing and CP violation. We do not have 
sufficient statistics to obtain useful measurements of the mixing and inter- 
ference CP violating parameters Am and sin0 (although they are left free 
in the fit). Thus, the primary results are the DCSD CP violating parameter 
Ad and the mixing parameters. The Fit D and D' results come from a fit to 
the combined and sample using Eq. 2. The global minimum for Fit D 
occurs at a highly unphysical x''^ = —0.52% and y' = 6.6% which makes the 
nominal branching ratio uninteresting. The minimum in the physical plane 
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M(D°) — WS GeV/c^ M(D°) — WS GeV/c- 




Fig. 2. Projection onto M(L>°) and (3(1?*) of both the fit (Fit C) and data for 
wrong-sign events. The left plots are for all events while the events in the right 
plots have cuts in the non-plotted variable. 

(Fit D') occurs at x'^ = 0.023% and y' = —2.6% with an increase in the 
log-likelihood fo 0.94 compared to Fit D. The branching ratio and CP asym- 
metry for the physical fit are found to be Rd = (0.521 Ig^gg ± 0.076)% and 
Ad = 0.13tQ;26 iO.lO. The statistical and statistical+systematic contours are 
determined in the same manner as in Fit C, as are the 95% CL limits on 
the mixing parameters. A mini Monte Carlo study using the best fit function 
to create fake data sets was employed to check the 95% CL contour. Only 
92.8 ± 0.8% of the fake data sets had fitted x',y' values inside the contour. 
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Fig. 3. Projection onto t{D^) of the fit (Fit C) and data for wrong-sign events. The 
left plot is for events in the M(D^) and Q{D*) signal region while the right plot 
shows the contributions to the wrong-sign signal. The first bin is low due to the 
necessary inefficiency of vertex separation requirements. The mixing component is 
negative due to the unphysical value of x'^ = —0.006 



Therefore the contour was increased to include 95% of the fake data sets by 
increasing A log£ from 3.00 to 3.50. The Fit D' contours are shown in Fig. 4. 
The improvement in — log£ of Fit D (D') compared to Fit B is 2.51 (1.57). 

The last fits, E and E', fit the and samples separately using Eq. 1. 
From these fits we can calculate Rd = \J RdRd = (0.518^o;|45 ±0.076)% and 
Ad = 0.12i[J:l ± 0.10. Contours at the 1 - VcTOS = 77.6% (Alog£ = 1.50) 
are obtained for each flavor. Every point in the contour is combined with 
every point in the D° contour. The 95% CL x', y' contour is defined by the 
outer edge of the resulting points. The systematic error contour is constructed 
in a fashion similar to Fits C and D'. The individual and contours are 
inflated to Alog£ = 1.50 + 0.12 = 1.62 and the combined contour found as 
before. A mini Monte Carlo study using fake data sets created from the best 
fit to the data was performed to check the 77.6% CL. The (-0°) contour 
was found to encompass significantly less (more) than 77.6% of the fake data 
sets. Therefore, the contours were each adjusted to include 77.6% of the fake 
data sets by changing Alog£ from the default 1.50. After this adjustment, 
the constructed 95% CL limit was found to be consistent with including 95% 
of the fake data sets. The Fit E' contours are shown in Fig. 4. Limits on the 
mixing parameters are set as in Fits C and D' by projecting the contour onto 
the axes. 
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Fig. 4. 95% CL x', y' contours for CP conserving and CP violating fits with x'^ 
constrained to be greater than zero. The inner curves show the 95% CL contours 
including statistical uncertainties only while the outer curves also include systematic 
uncertainties. 



6 Conclusions 



A compendium of the results from this analysis can be found in Tables 2 and 
3. The contours obtained in this analysis [9] agree better with recent results 
from BABAR [4] and BELLE [10] than with older results from CLEO [11]. 
The no mixing, no CP violation result agrees with the previous FOCUS re- 
sult obtained from the same parent data sample but utilizing a very different 
technique [2]. The two different techniques used to obtain the mixing with 
CP violation results, Fits D^'^ and E'^'\ return essentially the same answer. 
Fit D^'^ is to be preferred, however, since it retains all correlations in the fit. 
The results reported here represent the best published charm mixing limit 
from a fixed-target experiment. In some respects, this measurement is com- 
plementary to that obtained from the cleaner, higher statistics e^e~ collider 
experiments. FOCUS has the world's most accurate lifetime measurements for 
the six charm species it has measured. In addition, the FOCUS lifetime resolu- 
tion is far superior to that of the e+e^ collider experiments. These strengths, 
along with significant differences in production and reconstruction, result in a 
very different analysis from those of the e~^e~ collider experiments, providing 
a valuable check on those results. 
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Table 2 

Summary of results for Fits A-D, described in the text. 

Measurement Value 95% CL limit 

Fit A: No mixing, CP conserving 
Rws (0.4291°:°^? ±0.027)% 

Fit B: No mixing, CP violating 
Rws (0.429 ± 0.063 ± 0.028)% 

Ad 0.18 ±0.14 ±0.04 -O.IK yl^ < 0.48 

Fit C: Mixing, CP conserving, no x'^ constraint 
Rd (0.3811°:J^^± 0.092)% 

x'^ -0.059% 
y' 1.0% 

Fit C: Mixing, CP conserving, x'^ > constraint 
Rd (0.395 1 ± 0.069)% 

Rm < 0.63% 

x'^ 0.00% < 0.83% 

y' 0.5% -7.2% < y' < 4.1% 

Fit D: Mixing, CP violating, combined and fit, no x'^ constraint 
Rd (0.255lg:}^^± 0.132)% 

Ad 0.06lg:^^±0.13 
x'^ -0.52% 
y' 6.6% 

Fit D': Mixing, CP violating, combined and D° fit, x'^ > constraint 

Rd (0.517l[]:|5^± 0.076)% 

Ad 0.13tg|f ±0.10 -0.42 < < 0.81 

Rm < 0.63% 

x'^ 0.023% < 0.80% 

y' -2.6% -11.2% <y' < 6.7% 
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Table 3 

Summary of results for Fit E, described in the text. 

Measurement Result Result 

Fit E: Mixing, CP violating, separate and fits, no x'^ constraint 

R% (0.275lg;|^± 0.117)% (0.458 1 ± 0.099)% 

Rd (0.355l°:}^|± 0.107)% 

Ad -0.25 ±0.46 ±0.14 

x'^' -2.12% 0.61% 



y'^ 9.88% -4.43% 







Fit E': Mixing, CP violating, separate and fits, x'^ > constraint 

R^ (0.5861°;}^? ±0.079)% (0.458t°:|^ ± 0.099)% 

Rd (0.5181°;}^^ ^0.076)% 

Ad 0.12l°;^9±0.10 

Ad -0.46 <Ad< 0.72 @ 95% CL 

Rm < 0.61% @ 95% CL 

x'^' 0.00% 0.61% 

y'^ -1.02% -4.43% 

x'^ 0.023%, < 0.77% @ 95% CL 

y' -2.6%, -11.0% <y' < 6.6% @ 95% CL 
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